Exercise #8 — Descriptive Spatial
Statistics (Surfer)
Introduction to GIS Modeling, GEOG 3160,
University of Denver, Winter 2001
Name _____<enter your name>_____
Date _____<enter the date>______
From
the class website, download the data file Elev_pts_surfer.xls and place
it in
C:\Program Files\Golden Software\Surfer\Samples
Access
the Surfer demo system by double-clicking on C:\Program
Files\Surfer\Surfer.exe.
1) Derive an IDW surface of the elevation data.
a)
Generate a map surface…
ü From the main menu click Grid
à Data, then select the Elev_pts_surfer.xls
file in the …\Samples subdirectory.
ü In the Grid Line Geometry
box, change the # of lines to 25 in both the X and Y directions.
ü Click on the folder icon in
the Output Grid box, specify Elev_IDW1.grd as the output file
name, then press Save.
ü In the Gridding Method
box, specify Inverse Distance to a Power as the gridding method.
ü Click OK to generate
an IDW surface using the default specifications.
b)
Create a wireframe plot…
From
the main menu click MAP à Wireframe à OPEN à OK to create a wireframe plot
of the Elev_IDW1 surface.
c)
Create a contour map…
From
the main menu click MAP à Contour à New Contour map à OPEN à OK to create a wireframe plot
of the Elev_IDW1 surface.
In
the Contour Map Properties dialog box, check the Fill Contours
box on the options tab. Click on the Levels
tab, click on the Fill button, and click on the Foreground Color
button to pop-up the Color Spectrum dialog box.
Click on the left arrow of the spectrum and assign the color green; click on the right arrow and assign red. Ctrl/click in the middle of the spectrum and assign yellow. Click OK, OK and OK to generate the color filled contour plot.
d)
Graphically overlay the contour and wireframe maps…
Shift/click
to select both
the contour and the wireframe plots.
Then on the main menu, click Map
à Overlay Maps.
Resize and position the
colorful plot in the top third of the worksheet.
2) Repeat the “generate a map surface…” steps
substituting Elev_IDW2 for the output file name.
Press
the Options button on the right of the Gridding Methods box.
Options dialog box for IDW gridding
method
Set
Power to 4, Ratio 2, Angle 45 and press OK, OK
to generate a surface using different gridding method specifications.
Repeat
the “create a wireframe plot…, create a contour plot…, and graphical overlay…”
steps and resize/position the plot to the middle third of the worksheet.
3) Generate a third surface and resize/position
its plot in the bottom third of the worksheet that has the specifications: ELEV_IDW3,
Power 4, Ratio 2, Angle 45, Search Factors 1,1,1 and 1 (see
below).
Search Tab for gridding
Question 1.
Screen-grab your work sheet with the three plots and embed below…
Briefly describe the effects of changing the gridding method
specifications. Your discussion should
include explanation of IDW algorithm and the mathematical impact of changing
the power, anisotropy and search factors.
Referring to the screen grab extend the discussion to how the
changes you introduced affected the interpolated surfaces.
1) Create a new plot window by clicking on the New
button in the upper left corner of the main tool bar and selecting Plot
Document. Derive a Krig surface of
the elevation data as outlined below.
a)
Generate a map surface…
ü From Surfer’s main menu
click Grid à Data, then select
the Elev_pts_surfer.data file in the …\Samples subdirectory.
ü In the Grid Line Geometry
box, change the # of lines to 25 in both the X and Y directions.
ü Click on the folder icon in
the Output Grid box, specify Elev_Krig1.grd as the output file
name, then press Save.
ü In the Gridding Method
box, specify Kriging as the gridding method.
ü Click OK to generate
a Krig surface using the default specifications.
Click OK to the warning (any idea what it means?)
b)
Create a wireframe plot…
ü From the main menu click MAP
à Wireframe à OK to create a wireframe plot of the Elev_Krig1
surface.
ü Resize and position the
plot in the top third of the worksheet.
2) Repeat the “generate a map surface…” steps
substituting Elev_Krig2 for the output file name.
Press
the Options button on the right of the Gridding Methods box.
Options dialog box for Kriging
gridding method
Set
Kriging Type to Block, Drift Type to Quadratic and press OK,
OK to generate a surface using different gridding method specifications.
Repeat
the “create a wireframe plot…” steps for both of the new grids. Resize and position the interpolated
surface to the middle third of the worksheet.
3) Generate a difference surface between the
pair of Krig surfaces by Gridà Math and specifying
the Elev_Krig1.grd file as “Input Grid File A.” Specify Elev_Krig2.grd as input grid
file B. Specify Elev_%Krig1-2.grd
as output grid file C. Enter the grid
math function as C=( ( A – B ) / A ) * 100.
Plot,
resize and position the percent change surface (Elev_%Krig1-2.grd) to
the bottom third of the worksheet.
Question 2.
Screen-grab your work sheet with the three plots and embed below…
Make a general statement comparing the results of the IDW and
Kriging interpolations.
Which do you think is the “best” (what considerations lead to this
opinion)?
Briefly describe the effects of changing the Kriging gridding
method specifications. Your discussion
should include a general explanation of Kriging approach and the impact of
changing the type and drift options.
What is the data range of %change values?
Comment on the spatial pattern of the %change surface and suggest
any explanations for the differences.
From
Surfer’s main menu click Open à Elev_Pts_Surfer.xls
(Excel file), OK.
Highlight
the data in column C (Value) then select Data à Statistics to pop-up the statistics
dialog box.
Statistics dialog box
Click
OK to generate a standard report.
Click Close once you have reviewed the report.
Select
Data à Statistics and click
the check boxes for the additional statistics noted in the following table.
|
Sample Data
|
Mapped Data (IDW Surface) |
Number
of values |
|
|
Sum |
|
|
Minimum |
|
|
Maximum |
|
|
Range |
|
|
Mean |
|
|
Median |
|
|
First
quartile |
|
|
Third
quartile |
|
|
Standard
error of mean |
|
|
95%
confidence interval |
|
|
99%
confidence interval |
|
|
Variance |
|
|
Average
deviation |
|
|
Standard
deviation |
|
|
Coefficient
of variation |
|
|
Coefficient
of Skewness |
|
|
Coefficient
of Kurtosis |
|
|
Kolmogorov-Smirnov
stat |
|
|
Complete
the “Sample Data” statistics column in the above table using the
information in extended report.
From
the current window, switch to the Plot 2 window by clicking Windowà Plot 2. Now click Fileà Open à Elev_IDW1 (Surfer file), Open. Inspect the grid values by clicking on the
small crosses superimposed on the contour map.
Close the data inspection window (click on the small “X” in the upper
right corner of the plot).
From
Surfer’s main menu click Grid à Convert à Elev_IDW1.grd, Open to open
that file (map surface) for converting to a different data type.
Specify
Elev_IDW1_xyz as the file name and ASCIIXYZ (.dat) as the file
type.
BE SURE to the specify
ASCIIXYZ as the file type before pressing the Save button.
Click
Open à Elev_IDW1_xyz (.dat
file), OK. This will open the
grid file as a Worksheet.
Highlight
the data in column C, select Data à Statistics, and generate the descriptive statistics needed to
complete the left column of the above table.
To get the definition of each of the
statistics, press Helpà Contents then Search. Choose “Statistics Results” as the topic,
then click on the “Statistics” link as shown in the figure. Review the definition of each of the
descriptive statistics you entered into the table.
Embed below the definitions for any of the statistics that might be unfamiliar to you…
Question 3.
Comment on the comparison of the statistics for the Sample Data and the
Mapped Data.
Would you expect the statistics to be identical, nearly the same,
loosely related, or unrelated for most interpolated surfaces? Justify your answer.
Would you expect the Elev_Krig1 mapped data’s statistics to be
more, less or unpredictable in its comparison to the sample data’s
statistics? Justify your answer.